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We conjecture that thermalization following a quantum quench in a strongly correlated quantum 
system is closely connected to many-body derealization in the space of quasi-particles. This sce- 
nario is tested in the anisotropic Heisenberg spin chain with different types of integrability-breaking 
terms. We first quantify the deviations from integrability by analyzing the level spacing statistics 
and the inverse participation ratio of the system's eigenstates. We then focus on thermalization, by 
studying the dynamics after a sudden quench of the anisotropy parameter. Our numerical simula- 
tions clearly support the conjecture, as long as the integrability breaking term acts homogeneously 
on the quasiparticle space, in such a way as to induce ergodicity over all the relevant Hilbert space. 
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I. INTRODUCTION 

The understanding of ergodicity and thermalization in 
quantum systems is one of the most intriguing problems 
in quantum physics. Starting with the 1929 paper of 
John von Neumann 1 -, various attempts have been made 
towards the characterization of ergodic behavior in quan- 
tum systems^— and the establishment of a link with the 
notion of quantum chaos 4 ^. Theoretical interest in these 
issues resurfaced periodically^ until very recently, when 
an experimental study of the non-equilibrium dynam- 
ics of a quasi-one-dimensional condensate clearly demon- 
strated the lack of thermalization/ergodicity in a quan- 
tum many-body system 7 .. The attribution of this obser- 
vation to quantum integrability generated a lot of interest 
on its connections with ergodicity and thermalization in 
strongly-correlated quantum systems^. 

The simplest setting to study the relaxation of many- 
body systems is to consider an abrupt change in time 
of one of the control parameters, i.e., a quantum quench. 
At long times after the quench, the lack of thermalization 
in an intcgrable system can be seen as a consequence of 
the sensitivity to the specifics of the initial state encoded 
in the values of the constants of motion^. This lead to 
the proposal of describing the time-averaged steady state 
reached after a quench by keeping track of the initial 
value of all the constants of motion through a general- 
ized Gibbs ensemble^, whose conditions of applicability 
and drawbacks have been extensively tested^. In turn, if 
the system is far enough from the intcgrable limit, ther- 
malization is generally expected to occur, as numerically 
confirmed in many circumstances^—. 

This qualitative picture, although very appealing, 
leaves a number of important questions unanswered. It 
is not yet clear what is the nature of the integrable / non- 
integrable transition. Moreover, as it was shown in 
Ref. Il3l . it appears that even an integrable system could 
look "thermal", depending on the observable which is 
analyzed. Operators which are non-local in the quasi- 
particles of the system may behave thermally, while local 



operators do not. How to reconcile all these observations 
under a unifying framework? 

The purpose of this paper is to show that the under- 
lying mechanism governing the thermalization of many- 
body systems (and its relation to integrability) is that of 
many-body localization in Fock spaced— . We organize 
the paper as follows. In Sec. [HI after explaining in details 
the key concept of many-body localization, we discuss 
qualitatively our conjecture on thermalization following 
a quantum quench, i.e., on the role played by many-body 
localization. The model under investigation is then in- 
troduced in Sec. IIIII In the following sections we discuss 
our results: first of all we address the spectral properties 
of the model and use them to characterize the localized 
and dclocalizcd regimes (Sec. IIV[) ; we then focus on the 
quench dynamics, providing evidence of the connection 
between derealization and thermalization according to 
standard statistical mechanics predictions (Sec. |Vj) . Fi- 
nally, in Sec. I VII we draw our conclusions. 



II. MANY-BODY LOCALIZATION AND 
QUANTUM QUENCHES 

In this section we discuss a qualitative scenario con- 
necting the physics of thermalization after a quantum 
quench to the phenomenon of many-body localization. 
In order to do so, let us first lie down a few basic facts 
about many-body localization, as originally discussed in 
the context of transport of interacting electrons in ran- 
dom potential o 14 ' 16 . In absence of electron-electron in- 
teractions, the physics of disordered electron systems can 
be understood in terms of the standard Anderson local- 
ization phenomenology: extended wave functions corre- 
spond to finite zero-temperature conductivity while lo- 
calized states correspond to vanishing conductivity. The 
Anderson localization-delocalization transition is there- 
fore naturally associated to a metal-insulator quantum 
phase transition. Notice that since localized and ex- 
tended states cannot mix in the spectrum, the latter is a 
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sequence of bands of extended and localized states sepa- 
rated by mobility edges. 

The nature of the spectrum and of the eigenstates 
can change drastically if electron-electron interactions 
are taken into account. In particular, it has been re- 
cently showr>i^ that even when all single-particle states 
are localized, the presence of electron-electron interac- 
tions and inelastic collisions can result in an insulator- 
to-metal transition as the temperature of the system is 
raised above a certain critical value T c . Such a phase 
transition can be thought of as a many-body localization- 
delocalization transition, occurring at the level of many- 
body eigenstate a 14 ' 16 . Indeed, the presence of the many- 
body localization transition at finite temperature implies 
the existence of a many-body mobility edge at an energy 
scaling extensively with system size and separating local- 
ized many-body states, at low energies, from extended 
many-body states, at higher energies^. 

Many-body localization is a rather general concept 
which does not necessarily refer to real space. A stan- 
dard example in this sense comes from the physics of 
quasi-particle relaxation in quantum dots. Here the con- 
cepts of localization and derealization find their natural 
applicability in Fock space, where all many-body eigen- 
states are defined^. Similar ideas were also employed 
to analyze the mixing of vibrational modes due to an- 
harmonicity in molecules^. Along these lines, it was 
recently realized^— that the many-body localization- 
to-delocalization transition discussed above should be 
deeply connected to the main subject of this work, the 
physics of integrability-brcaking. More specifically, let 
us think of an integrablc model (having well defined 
quasi-particlcs) as a multidimensional lattice in which 
each point, identified by the occupations n(k) of the 
various quasi-particle modes, represents an eigenstate 
\^>a) = \{ n a{k)}) (see Fig. [T|). The space of these states 
(the quasi-particle space) is an obvious generalization of 
the standard Fock space. As long as states are localized in 
quasi-particle spaced, one expects the system to behave 
as integrable: any initial condition spreads into few sites, 
maintaining strong memory of the initial state. On the 
other hand, once a strong enough integrability-breaking 
perturbation hybridizing the various states \n a (k)) is ap- 
plied, a consequent derealization in quasi-particle space 
will occur (see Fig. [1]) . A tendency towards ergodicity is 
expected in this case. 

In this paper we aim at establishing a close connec- 
tion between the physics of the localization-dclocalization 
transition/crossover, occurring in quasi-particle space in 
the presence of an integrability-breaking term, and the 
physics of thermalization. In order to do so, we fo- 
cus on a specific class of non-equilibrium protocols on 
which thermalization can be studied, the so-called quan- 
tum quenches. In the present context they are defined 
through the time-dependent Hamiltonian 

n(t) = n [g(t)] + K lb , (i) 
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Figure 1: (color online). A cartoon of the quasi-particle 
space. For an integrable model all states, represented by 
the occupations of quasi-particles {n(k)}, are localized. An 
integrability-breaking perturbation introduces hopping ma- 
trix elements V among different sites, which hybridize, pro- 
vided \E{{n'(k)}-E({n"(k)}\ < V. For strong perturbations 
this may lead to derealization of wave functions among all 
points in quasiparticle space in a microcanonical energy shell. 

where: 

««>-{?£ :>»• < 2 » 

The time-dependent part of the Hamiltonian %o [<?(£)] is 
integrable, while Hib is the integrability-breaking term. 
We then ask ourselves the following question: what are 
the characteristics that many-body eigenstates should 
have in order for the system to thermalize and behave 
ergodically? 

The answer to an analogous question for semiclassical 
quantum chaotic systems was conjectured by M. Berry 
in 1977— and later employed by M. Srednicki to discuss 
thermalization in a (non-integrable) gas of interacting 
particles^. Inspired by these seminal papers, we propose 
that for generic many-body systems thermalization will 
occur whenever the eigenstates of the system become dif- 
fusive in microcanonical shells defined in quasi-particle 
space. These diffusive states correspond to the intuitive 
expectation that in an ergodic state any initial state is al- 
lowed to diffuse into all states in a micro-canonical energy 
shell, generating a cascade of all possible lower energy ex- 
citations^. The purpose of the remaining sections is to 
test this proposal on the dynamics of a concrete inte- 
grable model. 

III. THE MODELS 

In order to corroborate the scenario proposed above, 
we now study in detail the dynamics after a quench of 
an anisotropic Heisenberg chain subject to various forms 
of integrability-breaking perturbations. The anisotropic 
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Heisenberg (XXZ) spin-1/2 chain is denned by: 
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where of (a = x,y,z) denote the spin-1/2 Pauli ma- 
trices on site i, J the coupling strength, while J z the 
z-anisotropy intensity. Units of H = ks = 1 are used 
throughout the paper; J = 1 is taken as the energy scale. 
This Hamiltonian is integrablc by Bethc Ansatz, and ex- 
hibits two gapped phases, ferromagnetic (J z < —1) and 
antifcrromagnetic ( J z > 1), separated by a critical region 
— 1 !_ Jz < 1; with J z -dependent critical exponents^ and 
quasi-long-range-order in the xy spin-plane. 

As for the integrability-breaking perturbation, we con- 
sider different cases, with or without disorder terms in 
the Hamiltonian, which can be expressed in the form: 



(4) 



where is the amplitude (possibly site-dependent) of 
an additional few-body term 0j. This few-body term 
may act on a single site i (e.g., onsite magnetic field), or 
on a few sites centered around i (e.g., nearest or next- 
to- nearest neighbor couplings). In particular, we break 
integrability by either adding: (/) a random magnetic 
field in the z-direction; (//) random J z couplings; (///) 
random or (IV) uniform next-nearest neighbor zz cou- 
plings, according to: 
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For cases (I)- (I I)- (1 1 1) integrability-breaking is induced 
by the disorder, hi 6 [—1,1] being random numbers, 
while in case (/V) disorder is not invoked. In general, 
disordered systems allow for a better statistical analy- 
sis, due to the possibility of averaging over randomness. 
One might argue that such averages are strictly required 
in order to reproduce our findings about thcrmalization. 
This is not the case, since we have found analogous qual- 
itative conclusions in all the four cases discussed above: 
integrability-breaking is the only crucial requirement for 
our mechanism of thcrmalization to set in. 

We first address the spectral properties of the model, 
and subsequently consider a sudden quench of the 
anisotropy parameter J z = g. The total magnetization 
S z = Ei a l is a conserved quantity, hence we restrict to 
the sector S z = 0. Nonetheless, due to the involvement 
in such non-equilibrium dynamics of a considerable part 
of the spectrum, standard (both analytic and numerical) 
rcnormalization group techniques are eventually doomed 
to failure. We therefore resort to exact numerical diago- 
nalization of systems with up to 16 spins. 



While the zero-temperature phase-diagram in presence 
of disorder is well established [25|, the high-temperature 
phase-diagram has been conjectured to be composed of 
two phases, a non-ergodic many-body localized phase 
(in real space) at A > A crlt , and an ergodic one at 
A < A crit ; in the case (I), A crit - 6 ~ 8 at J z = 1 
(in our units ) 18 ' 19 . The results presented below indicate 
the presence of a second non-ergodic localized phase (in 
quasi-particle space) for A close to zero that crosses over 
to the ergodic phase upon increasing A. The fate of this 
crossover in the thermodynamic limit and the eventual 
value of the critical A* are yet to be determined^. 

We start by characterizing deviations from integrabil- 
ity in terms of the many-body level statistics and of 
the properties of the eigenstatcs. A well defined tran- 
sition from Poisson (Integrablc) to Wigncr-Dyson statis- 
tics (non-Intcgrablc) is closely associated to the local- 
ized/diffusive character of eigenstates in quasi-particle 
space. Using this characterization, we then show that the 
non-thermal-to-thermal transition in the dynamics is di- 
rectly connected to the localization/delocalization tran- 
sition in quasi-particle space. In particular, by looking 
at the asymptotics of spin-spin correlation functions, we 
discuss how thcrmalization is linked to the emergence of 
diffusive eigenstates in quasi-particle space. This also al- 
lows us to discuss, in a broad context, the relationship be- 
tween locality of observables in quasi-particle space and 
the corresponding behavior. 



IV. SPECTRAL PROPERTIES 

Let us first concentrate on the spectral properties of 
the Hamiltonian for a given value of the anisotropy g = 
J z . In the following we show data for J z = 0.5. We have 
explicitly checked that, changing J z to a different value 
within the critical region of the XXZ model (\J Z \ < 1), 
does not qualitatively affect the scenario discussed and 
our conclusions. 



A. Level Spacing Statistics 

The statistics of the energy levels represents a key fea- 
ture of the spectrum of a generic quantum system, since it 
is a good indicator of the presence of integrability. Both 
in semiclassical and in many-body systems, integrable 
systems have levels that tend to cluster, eventually cross- 
ing when a parameter in the Hamiltonian is varied. On 
the other hand, in non-intcgrablc systems the levels are 
correlated in such a way as to avoid crossings. A quan- 
titative way to characterize these tendencies is through 
the Level Spacing Statistics (LSS)2l, i.e., the probability 
distribution P(s) that the energy difference between two 
adjacent levels s n = E n+ \ — E n (normalized to the av- 
erage level spacing) falls in the interval [s, s + ds]. In a 
typical integrable system one finds a Poissonian (P) LSS: 

P P (s) = e- s . (5) 
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On the other hand, for non-integrable systems one ex- 
pects random matrix theory to apply, leading to a 
Wigncr-Dyson (WD) distribution, where level repulsion 
shows up in lim s _j.o P(s) ~ s 7 . More specifically, for 
systems as the one considered here, which preserve one 
anti- unitary symmetry (invariance under time-reversal), 
the statistics is given by a Gaussian Orthogonal Ensem- 
ble^ (at low energy spacings one has the characteristic 
behavior 7=1): 



Pvm{s) = -y e 



(6) 



In our case the system undergoes a transition from 
Poissonian [Eq. (0] to Wigner-Dyson [Eq. flB])] LSS upon 
increasing the non-integrable perturbation A, which for 
finite-size systems takes the form of a smooth crossover. 
This can be faithfully quantified by means of the Level 
Spacing Indicator (LSI) 77: 



n 



Io°[Pwd(s) 



P P (s)}ds ' 



(7) 



where P(s) is the probability distribution function of 
the level spacing between neighboring levels, while sq ~ 
0.4729 is the first intersection point of Pp(s) and Avd(s). 
The LSI is zero for systems with a Poisson distribution 
Pp of the spacings and one if the distribution is Wigner- 
Dyson Pwd- Below we will study the LSI in two different 
ways: (i) as a function of the energy eigenvalues in indi- 
vidual microcanonical shells [ry w (£?)], using Eq. (0 with 
P\e,e+w] ( s ) as th e level statistics computed in the en- 
ergy window [E 7 E + W}; (ii) cumulatively \r] c (E)], as a 
function of the energy eigenvalues below a given thresh- 
old, with P[e ,e]{s) as the level statistics of eigenvalues 
with excitation energy less than E, with respect to the 
ground state energy Eq. 



1. Results for the LSI 

The LSS in the XXZ model with an integrability- 
breaking perturbation has been the subject of various 
studies in the literatur o 28 i 29 . Here we are not interested 
in a complete characterization of it, but rather on elu- 
cidating under which conditions, and in which regions 
of the energy spectrum, the model behaves according to 
the WD statistics, i.e., rj « 1 according to our defini- 
tion. We point out that the WD distribution of Eq. © 
is obtained for non-integrable systems with only a time- 
reversal symmetry. In all our simulations we considered 
open boundary conditions, fixed the sector of zero mag- 
netization, and added a very small magnetic field on the 
first site of the chain, in such a way as to work in a sub- 
space without any other unwanted symmetry. 

All the four types of intcgrability-breaking perturba- 
tion behave quite in the same way, the only difference 
being for case (/V), where fluctuations are more consis- 
tent, due to the absence of disorder averaging (to reduce 
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Figure 2: (color online). Level spacing indicator for the XXZ 
model © with non-integrable perturbations (I) (left panels) 
and (II) (right panels). The upper frames show -q„(E), while 
the lower ones display rj c (E). Data are for L — 14 sites with 
different values of the integrability-breaking perturbation A. 
The LSI ?7w (E) is evaluated in a microcanonical shell of width 
W = 2. Averages are performed over 10 3 disorder instances. 
In order to exactly recover Poisson and GOE statistics in the 
two integrable and non-integrable limits, we performed an 
unfolding of the energy spectrum for each instance, according 
to standard techniques adopted in quantum chaos^. 



fluctuations, one should consider energy spectra of larger 
systems; however, the exact diagonalization technique in- 
trinsically imposes severe size limitations) . As an explica- 
tive example, in Fig. [5] we plot both rj w (E) and rj c (E) for 
the XXZ model with a random z-field, (J) (left panels), 
and with random J z couplings, (II) (right panels). We 
observe that, fixing the system size, if A is progressively 
increased, the value of 77 also increases, until it reaches, 
in the middle of the energy band, a value close to 1 (for 
A ~ 1, in our units and at L = 14). For A > 1, 77 de- 
creases again towards small values, since for A> J z the 
system turns into a trivial classically integrable model 2 ^. 
Only in the middle of the spectrum the system appears 
to exhibit level repulsion, while this is not the case in the 
low- or high-energy spectrum^. This is more evident 
from the cumulative LSI rj c ; here one can notice that, for 
sufficiently strong perturbations and at low energies, r] c 
is an increasing function of E, until it saturates around 
its maximal value. 



B. Inverse participation ratio 

After a spectral characterization of integrability 
through the LSS, we come back to the characterization of 
the eigenstates. In the region where a Poissonian LSS is 
observed, the eigenstates are expected to be localized in 
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quasiparticlcs space, while if WD is seen, the eigenstates 
are expected to be delocalized. The proper tool to quan- 
titatively characterize the properties of the eigenstates 
and their derealization is the so called Inverse Partici- 
pation Ratio (IPR)2L2a. The IPR on a normalized pure 
state \ip) is a basis-dependent quantity, defined by: 



N 



am 



(8) 



where {\n)} is the reference basis of the Hilbert space. If 
a state is a uniform superposition of n st basis states, the 
corresponding contribution to £ is of order n s t- 

We will focus on the IPR of the system eigenstates, 
evaluated of two types of basis: (i) the site (S) basis 
\ n s) = \o~i ■ ■ ■ (Tl) = ±1), composed by the eigen- 
states of of, which is often referred to as the "compu- 
tational basis"; (ii) the integrable (I) basis, composed 
by the eigenstates of the integrable model <j3j) in absence 
of the perturbation terms: A = 0. Analogously to the 
LSI, we can compute the IPR over microcanonical shells 
around a given energy value E. Notice that, if an eigen- 
state is localized in quasi-particle space, we expect the in- 
verse participation ratio computed in the integrable basis 
to be £r — 0(1). Conversely, if an eigenstate is a diffu- 
sive superposition with random phases and similar am- 
plitudes of N eigenstates of the integrable model, then 
£j ~ N , Below we will use these facts to characterize 
localization and derealization in quasi-particle space. 



1. Results for the IPR 

In Fig. |3] we show the inverse participation ratio for 
the XXZ model with the integrability-brcaking terms (J) 
and (II) (as for the LSI, results are not qualitatively dif- 
ferent if different perturbations are considered) . Looking 
at the IPR in the site basis (lower panels) , as long as A is 
increased we observe a general tendency to a localization 
(the IPR peak value decreases) . This is coherent with the 
fact that the states of the computational basis arc exactly 
the eigenstates of the system for A ^> J z . On the other 
hand, as depicted in the upper panels, the IPR in the 
integrable basis behaves rather differently. In particular, 
it provides a clear signature of the fact that eigenstates 
delocalize with increasing values of the disorder A. 

In order to better understand the nature of the de- 
localization induced by the integrability-breaking term, 
let us compare £ with the number of states N\e e+W] 
in the relevant microcanonical shell [E,E + W], where 
W ~ V is the typical matrix element of the integrability- 
breaking perturbation (V « 2 A in this case). For small 
A, £/ -c N[ E ,e+w] (see Fig. upper panels), indicating 
that the eigenstates are still close to those of the inte- 
grable system and the degree of derealization of the sys- 
tem is very low. On the contrary, when A ~ 1 (Fig. 21 
lower panels) we observe that £/ ~ N^ E ^ E+W ^ . In this 
case, the perturbation is able to hybridize nearly all the 
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Figure 3: (color online). Inverse participation ratio for the 
XXZ model ((3J with non-integrable perturbations (/) (left 
panels) and (II) (right panels). Data are for L = 14 sites and 
different values of the perturbation strength A. The IPR is 
evaluated in a microcanonical shell of width W = 2A, in the 
integrable (upper frames) and in the site (lower frames) basis. 
Data are averaged over all the system eigenstates in the ap- 
propriate energy window, and over 10 2 disorder realizations. 



quasiparticles states within the microcanonical energy 
shell. As we will see in the next section, this is the key in- 
gredient for the system to thcrmalize. Notice that in this 
context the low-lying eigenstates are rather peculiar: this 
part of the spectrum, which contains very few states as 
compared to the center, has closely Poissonian statistics 
and is characterized by large fluctuations of statistical 
quantities. 



V. DYNAMICS AFTER THE QUENCH 

It is now time to use the information we obtained 
about spectral statistics and eigenstates to study the re- 
lation between many-body localization and thcrmaliza- 
tion. We will look at the dynamics following a sud- 
den quench of the anisotropy parameter g = J z from 
Jzo at t < to J z ^ Jjo at t > 0, as described in 
Eq. ([2]). As we did for the spectral properties, we will 
show data for systems where the anisotropy is quenched 
toward J z = 0.5. Different values of such J z do not 
qualitatively affect the scenario. The system is initially 
prepared in the ground state l^o) of PL(Jm), so that its 
(conserved) energy with respect to the final Hamiltonian 
K(J Z ) is E Q = (ipo\'H(Jz)\ipo)- For growing values of J z0 , 
the state IV'o) tends towards the classical antiferromag- 
netic Neel state, and Eq/L saturates to a constant value, 
slightly below the middle of the spectral band, thus im- 
plying that a quench generally involves only a fraction of 
the eigenstates of the hnal Hamiltonian. 



Figure 4: (color online). IPR in the integrable basis £7 at 
A = 0.1 (upper panels), and at A = 1 (lower panels), com- 
pared to the number of states N in an energy window of width 
W — 2A. Results are shown for the integrability-breaking 
terms (/) (left panels) and (II) (right panels). Notice how- 
ever that all other cases display an identical qualitative be- 
havior. In particular, the disorder is not required to observe 
the hybridization of the quasiparticle states within the micro- 
canonical energy shell for the case of large A. 

A. Effective temperature 

Contrary to local quenches, the work done on the sys- 
tem by changing the anisotropy from J z q to J z is exten- 
sive. It is then interesting to ask, after a quench involving 
an extensive injection of energy Eq — E gs oc L [E SB being 
the ground state energy of H(J Z )] , if the subsequent long- 
time evolution of the system is effectively described by 
an equilibrium dynamics governed by T-L(J Z ). In view of 
a plausible equivalence between a microcanonical (fixed 
Eq) and a canonical equilibrium description of such a 
long-time dynamics, it is meaningful to define, as in pre- 
vious instances^, an effective temperature T e g for the 
system out of equilibrium. We compute T e g by equating 
the micro-canonical energy Eq = (ipo\ / H(J z )\ipo) to the 
canonical ensemble average 

E = {U(J z ))t o{{ = Tr [p(T eS ) %{J Z )} , (9) 

where p(T ff) is the equilibrium density matrix at tem- 
perature T e ff. 

e -H{J*)/T eS 

= Tr[e ^^(.7,)/T cl , ] ■ ( 10 ) 

This temperature is eventually averaged over disorder re- 
alizations in the cases (I)-(I 1 1) . 

In Fig. [5] we show the effective temperature as a func- 
tion of the initial value of the anisotropy for a system of 
L = 12 sites, quenched toward J z = 0.5. As it is appar- 
ent, T e ff is monotonically increasing with \ J Z — J z o\- In 



Figure 5: (color online). Effective temperature in the XXZ 
model with an integrability-breaking perturbation, after a 
quench in the anisotropy parameter toward a value of J z = 
0.5. The values of J z o in the x-axis denote the initial 
anisotropies, and stand for different initial conditions [that 
is, the ground states of the Hamiltonian H(J z o)]. Data are 
for L — 12 sites; in all panels, except the lower right one, 
averages are performed over 200 disorder instances. 



the first three cases, the effective temperature saturates 
for large values of J z q, because the initial ground state 
|V'o) tends toward the antiferromagnetic Neel state (for 
J z ^> 1 and A < 1 the effective temperature is around 
T c s ~ 5, thus meaning that the states probed are located 
in the lower central part of the band)2I. 

B. Thermalization of correlation functions 

We are now ready to test the relation between delo- 
calization and thermalization after a quench by studying 
the long-time asymptotics of two-spin correlators, con- 
structed as expectation values of 

1 L 

< = 1 E ^- l)k,L ofaf , (a = x, z). (11) 
3,1=1 

In particular, we compare the expectation value in the 
canonical ensemble at the corresponding T e s'- 

(k) ee «> Toff = Tr [p(T eS ) <] , (12) 

with the asymptotic value that is reached after the 
quench, calculated from the diagonal ensemble^: 

n a Q (k) ee lim (mK\m) - E \d\ 2 (<t>iK\4>i) , (13) 

i 

where \ip(t)} = e~ lW ( Ji '*|?/>o} is the state of the system 
at time t, while c,; = (4>i\ipo) is the scalar product be- 
tween the state |^o) and the eigenstatcs \<f>i) of the final 
Hamiltonian %{J Z ). 
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The obscrvablcs we consider here correspond to two 
completely different scenarios in terms of the system 
quasiparticlcs, being nf a local operator while n| a non- 
local one. We recall that, in this context, local and non- 
local operators refer to the structure of their matrix el- 
ements on the basis of quasi-particles: local means that 
the operator couples a finite number of states, while non- 
local that it couples all states^. While correlators in the 
x-direction are always well reproduced by an effective 
thermal ensemble, correlators in the z-direction appear 
to be more sensitive to the breaking of integrability. This 
is seen quite clearly in Fig. [SJ where we plot the correla- 
tors n£ averaged in the diagonal (black circles) and in the 
canonical (red squares) ensembles, both along the x-axis 
(upper panels) and along the z-axis (lower panels). The 
parameters are chosen in such a way as to have the sys- 
tem close to integrability, with a significant derealization 
in Fock space still not present. Two different behaviors 
for n x (k) and n z (k) are apparent, with differences clearly 
emerging at the peaks k = tt, where boundary effects are 
less pronounced. One can qualitatively see that, while 
discrepancies between the two ensembles are well visible 
in n z (k), they are suppressed in n x (k). Therefore, in a 
quasi-integrable regime only n x behaves thermally, while 
n z does not. As stated above, this reflects the intrinsic 
difference between nonlocal/local operators with respect 
to the quasiparticles, which emerges only for the cases in 
which the system itself is not able to properly hybridize 
the quasiparticle states within the microcanonical shell. 
Here we stress however that the classification of the op- 
erators is in general a subtle issue. In the model we con- 
sidered in this work, it has been possible by analyzing 
the XX limit (J z = 0) and the low-energy sector of the 
critical phased. 

A quantitative measure of the degree of thcrmalization 
is given by the absolute discrepancy between the diagonal 
and the canonical ensemble predictions: 



5n a k = \n«(k) - nlJk)\ 



(14) 



In order to elucidate the drastically different behavior 
between integrablc and non-integrable systems, in Fig. [7J 
we plot Sn k at the peak k = 7r where discrepancies are 
larger, as a function of the disorder amplitude A. We 
observe that Sn x is more than one order of magnitude 
smaller than Sn^, indicating a closely thermal behavior, 
while 5n^ shows a sharp decrease as integrability is pro- 
gressively broken by increasing A, towards a minimum 
value at A. The scaling with the dimension L of the 
chain, shown in the upper left panel of Fig. ([7]), confirms 
our predictions. While the behavior of Sn x is indepen- 
dent on the system size, the decrease of <$n* as a function 
of A is more pronounced on increasing L. Due to the nu- 
merical limitations of exact diagonalization, we cannot 
rule out the possibility that, in the thermodynamic limit, 
the integrable to non-integrable transition for low pertur- 
bation strengths occurs at A* = in all the cases ana- 
lyzed here (I)-(IV). However, for the considered sizes, 
we found that the local observable 5n z thcrmalizcs the 




Figure 6: (color online). Comparison between the diagonal 
and canonical expectation value of the two-spin correlation 
function n x (k) (upper panels) and n z (k) (lower panels) as a 
function of the momentum k. Data are for a quench from 
JzO = 10 to J z = 0.5 and two kinds of integrability-breaking 
perturbation [(7) in left panels, and (77) in right panels], with 
disorder intensity A = 0.4. 
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Figure 7: (color online). Discrepancies Sn% (black circles) 
and 5ni (red squares) between the diagonal and the canonical 
ensemble predictions, for a quench from J z q = 10 to J z = 0.5. 
Data in the upper left panel (7) are for different system sizes 
as depicted in the caption, while in the other panels are for 
L — 12. In the cases of perturbations involving disorder, (7)- 
(77)- (777), averages over 200 instances are performed. 



best at around A ~ 1 for model (7). It is now crucial 
to observe that around this point, as noticed in the pre- 
vious sections, the diffusive nature of the eigenstates is 
pronounced. The upper right and the lower left panels 
seem to locate the optimal thcrmalization point for mod- 
els (II) and (777) at a slightly larger value of A, while 
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the lower right panel shifts it to slightly smaller values 
for model (IV). In all cases it is however true that at 
these points £ ~ N[e,e+w] i making a direct connection 
between good thcrmalization and diffusive nature of the 
eigenstates. Notice also that, for A > A, Sn^. neces- 
sarily has to increase again, in agreement with the fact 
that A ~ A is the point where the non-integrable be- 
havior is most pronounced and that for large values of A 
the system tends toward another integrablc limit^. In 
analogy with previous studies, the different sensitivity to 
intcgrability of correlators in different spin directions can 
be qualitatively understood as a consequence of the fact 
that <j z is a local operator in quasi-particle space while 
a x is a non-local one^. 



system thermalizcs, one should study the correspond- 
ing many-body localization/delocalization transition in 
quasi-particle space. Thcrmalization should occur when 
the relevant typical states spread diffusively on an expo- 
nential number of states lying in the microcanonical en- 
ergy shell. We point out that our picture is valid as long 
as the integrability-breaking term acts homogeneously in 
the quasiparticle space, in such a way as to induce er- 
godicity over all the relevant Hilbert space. For generic 
dynamic systems there may be regions of the phase space 
which are non chaotic, so that their quantum versions 
produce entropy at a non-uniform rate given by the lo- 
cal Lyapunov exponents^. In this case more complex 
scenarios for the approach to equilibrium may arise. 



VI. SUMMARY 

In conclusion, we discussed thermalization and 
integrability-breaking in the dynamics after a quench of 
a quantum XXZ Hcisenberg spin chain in presence of an 
integrability-breaking term. We have shown that, if one 
wants to know when and how an interacting many-body 
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